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Abstract. We present three-dimensional numerical simulations 
of the atmospheric dynamics of close-orbiting planets such as 
HD209458b. Our simulations show that winds of several kmsec"^ 
and day-night temperature differences reaching 500-1000 K are 
possible at and above the photosphere. The circulation takes the 
form of a broad superrotating (eastward) equatorial jet. At ~0.1- 
1 bar, the superrotation blows the hottest regions of the atmosphere 
downwind by ~ 60° of longitude, but at lower pressures the temper- 
ature pattern tracks the stellar illumination. We predict factors of 
several variation in the infrared flux received at Earth throughout 
an orbital cycle; if the photosphere is deep enough (>50-100 mbar 
pressure), the peak infrared emission should lead the time of sec- 
ondary eclipse by 10 hours or more. Dynamics plays a key role in 
shaping the spectrum, clouds, chemistry, and long-term planetary 
evolution. 

The past few years have witnessed many observations constrain- 
ing the physical properties of extrasolar giant planets with orbital radii 
less than 0.1 AU (the "Pegasi planets" or hot Jupiters). Eight such 
planets have been discovered to undergo transits. Two such planets, 
HD209458b and TrES-1, have also been detected in thermal emission 
during the secondary eclipse, and several useful upper limits on compo- 
sition, albedo, and thermal emission at various wavelengths have been 
achieved. This trend of detections is likely to continue. 

A knowledge of atmospheric dynamics will be crucial for under- 
standing these new observations. First, the interaction of dynamics 
with radiative transfer controls the temperature structure, which shapes 
the infrared spectrum and lightcurve. Most current radiative-transfer 
models adopt radiative-equilibrium conditions and make arbitrary as- 
sumptions about whether the absorbed stellar flux gets redistributed 



1 



2 



Showman and Cooper 



across the planet or heats only the dayside; however, dynamics can push 
the atmosphere far from radiative equilibrium, and the extent of heat 
redistribution must be calculated explicitly (and may depend strongly 
on height). Second, whether clouds exist depends on the temperatures 
and locations of ascent /descent, which is again controlled by the cir- 
culation. Cloudiness in turn determines the albedo, visible lightcurves, 
and — if high-altitude clouds form — causes masking of spectral lines. 
Third, circulation may lead to disequilibrium between CO and CH4, re- 
move condensable species (Na2S, CaTiOa) via nightsidc cold trapping, 
and cause other chemical effects. Fourth, the atmospheric heat engine 
produces enormous kinetic energy, which, if transported deep enough, 
may affect the interior evolution. This has been suggested as a possible 
mechanism for producing the large radius of HD209458b, for example 
(Guillot and Showman 2002, Showman and Guillot 2002). 

The intense starlight incident upon the surface of Pegasi planets 
leads to a deep radiative zone extending from the top of the atmosphere 
to pressures of ~ 1000 bars (e.g., Guillot ct al. 1996, Cuillot and Show- 
man 2002, Burrows et al. 2003, Chabrier et al. 2004), and any observ- 
able weather occurs in this radiative zone. The fast spindown times for 
Pegasi planets implies that these planets should be in near-synchronous 
rotation (3.5 days for HD209458b) (Guillot et al. 1996, Showman and 
Guillot 2002) . This rotation rate implies that rotation is important but 
not dominating: for kmsec"^ winds, the Rossby number is ~ 1. These 
estimates imply dynamical length scales (the Rossby deformation radius 
and Rhines length) of order a planetary radius. As a result, any jets and 
gyres that exist should be global in scale. This contrasts with the case 
of Jupiter, where these length scales are only 2-10% of the planetary 
radius and — as a result — the dominant jets and vortices are much 
smaller than a planetary radius. Pegasi planets should therefore have 
physical appearances that differ greatly from Jupiter and Saturn. 

Here we describe our recent work on the atmospheric circulation 
of Pegasi planets; the presentation describes and extends Cooper and 
Showman (2005), to which we refer the reader for details. We per- 
formed global, three-dimensional numerical fluid simulations using the 
ARIES/GEOS dynamical core (Suarez and Takacs 1995) in a domain 
extending from Imbar — 3kbar. The simulations adopted the primitive 
equations, which are a simplified form of the Navier-Stokes equations 
valid for statically stable atmospheres that are vertically thin compared 
to their horizontal extent. We used parameters for HD209458b and as- 
sumed that the planetary interior is in synchronous rotation with the 
orbital period. The nominal resolution is 72 x 45 in longitude and 
latitude with 40 vertical levels. In the simulations, the dynamics are 
driven solely by the imposed day-night heating contrast. Rather than 
solving the radiative transfer explicity, we adopted a thermodynamic 
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Figure 1.: Temperature (greyscale) and winds (arrows) at pressures of 2 
mbar, 220 mhar, and 20 bars at 5000 Earth days of simulated time. Peak 
winds are 9.2, O'^d 2.8kmsec~^ from top to bottom, respectively. 
Heating occurs on the dayside (longitudes —90° to 90°j and cooling oc- 
curs on the nightside (longitudes -180° to -90° anc? 90° to 180°J. The 
substellar point is at 0° latitude, 0° longitude. 
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heating rate (in Ksec~^) of (Teq — T)/rrad, where Tgq is the specified 
radiative-equilibrium temperature profile (hot on the dayside and cold 
on the nightside), T is the actual temperature, and Trad is the radiative- 
equilibrium timescale (a function of pressure) . The vertical structure of 
Teq and Trad were taken from Iro et al. (2005); the day-night difference 
in Teq was a free parameter that we varied from from 100 — 1000 K. 

Figure 1 shows the temperature (greyscale) and winds (arrows) 
for three layers (2 mbar, 200 mbar, and 20 bars from top to bottom, 
respectively) after a simulated time of 5000 Earth days. The imposed 
heating contrast leads to winds exceeding several km see" ^. By 5000 
days the simulation has approximately reached a statistical steady state 
at pressures less than 3 bars, although the winds continue to increase 
at deeper levels. At the top (2 mbar), the radiative time constant 
is ~ Ihour (much less than the advection time), so the hot regions 
remain confined to the dayside. The temperatures are in near-radiative- 
equilibrium, with day-night temperature differences of ~ 1000 K. At 200 
mbar, close to the expected photosphere if the planet lacks high-altitude 
clouds, a broad ~ 4kmscc~^ eastward supcrrotation develops. Here, 
the radiative time constant, ~ 10^ sec, is comparable to the time needed 
to advect air across a planetary radius. The circulation therefore blows 
the hottest regions of the atmosphere downwind from the substellar 
point by ~ 60° degrees of longitude. Temperature differences reach 
~ 500 K at this level. 
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Figure 3.: Infrared brightness of HD209458b as viewed at Earth during 
(a) transit, (b) one-quarter period after tranit, (c) secondary eclipse, 
and (d) one-quarter period after secondary eclipse. The planetary ro- 
tation axes are vertical, with the superrotating jet seen in Fig. 1 going 
from left to right in each panel. The smaller schematic globes in the 
top row illustrate the illumination of the planet by sunlight, as viewed 
at Earth, during these same phases. Time increases from left to right. 



The development of a broad, superrotating (eastward) jet, with 
large day-night temperature differences at p < 1 bar, is a robust fea- 
ture in all our simulations. We performed a variety of simulations 
using radiative-equilibrium temperature profiles from Burrows et al. 
(2003), Chabrier et al. (2004), or Iro et al. (2005); and using radiative- 
equilibrium day-night temperature differences of 1000, 750, 500, 250, or 
100 K. We even performed simulations whose initial condition contained 
a broad westward equatorial jet extending from ~ 2 bars to the top of 
the domain. All of these simulations developed strong eastward jets re- 
sembling that in Fig. 1. Fvirthcrmore, these results agree with Showman 
and Guillot (2002), who also obtained broad eastward jets in every one 
of their simulations using a different numerical code. This gives us con- 
fidence that eastward flow is a robust result, at least within the context 
of our adopted input parameters. What all these simulations have in 
common are short radiative time constants at pressures < 1 bar, which 
allow the development of longitudinal temperature variations that are 
essentially a large-amplitude thermal tide. Wc speculate that, as has 
been suggested for Venus, this tide induces the superrotation by pump- 
ing eddy energy and eastward momentum upward and equatorward. (In 
absence of such eddy effects, the equatorial flow would be westward.) 

Our results differ from the one-layer shallow-water calculations of 
Cho et al. (2003), which produce westward equatorial flow. However, 
shallow-water turbulence invariably produces westward flow even for 
planets such as Jupiter and Saturn whose equatorial jets are eastward 
(e.g., Cho and Polvani 1996, lacono et al. 1999; see Vasavada and 
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Figure 4.: Predicted infrared lightcurve for HD209458b assuming black- 
body emission from the 220 mbar level (left) or the 10 mbar level (right). 



Showman 2005 for a review). This feature seems to result from the 
exclusion of three-dimensional processes in the shallow- water equations. 

The patterns in Fig. 1 have implications for the infrared lightcurve 
of the planet throughout its orbit, as shown in Figs. 3 and 4. The globes 
in Fig. 3 (bottom row) show the infrared brightness at four phases as- 
suming the planet emits as a blackbody from the 220 mbar level. The 
globes in the top row show the illumination as viewed from Earth dur- 
ing these same phases. The key point is that, in the absence of winds, 
the temperature pattern would follow the illumination (i.e., the infrared 
appearance would also correspond to the top row of globes). The dif- 
ferences between the idealized illumination patterns and the simulated 
brightness patterns result solely from atmospheric dynamics. 

Figure 4 (left) shows the corresponding lightcurve assuming the 
photosphere is at the 220 mbar level, as might be expected for a cloud- 
free planet. Because the hot regions become offset from the substellar 
point at this pressure, the model predicts that the planet will radiate 
its maximum infared flux toward Earth ~ 14 hours before the secondary 
eclipse (rather than immediately around the time of secondary eclipse 
as would occur without winds) . This effect could allow an observational 
determination of wind direction — if the winds are eastward, the peak 
fluxes would lead the eclipse, whereas if the winds are westward, peak 
fluxes would lag the eclipse. 

Nevertheless, uncertainties exist regarding the emission level. The 
lower-than-expected Na (Charbonneau et al. 2002) and null detection of 
CO (Deming et al. 2005, Richardson et al. 2003) on HD209458b could 
result from a high-altitude cloud at a few mbar pressure. If the cloud 
optical depth exceeds unity, then the primary infrared emission to space 
occurs from the cloud altitude rather than from the deeper levels ex- 
pected for a cloud-free planet. Figure 4 (right) shows the lightcurve for 
the case of emission from the 10-mbar level. Because the radiative time 
constant is short at these pressures (Iro et al. 2005), the offset shown 
in Fig. 4 (left) has largely disappeared. The magnitude of the flux dif- 
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Figure 5.: Solid and dashed curves show the temperature profile at 
the suhstellar and antistellar points, respectively, for the simulation in 
Fig. 1. The dotted line shows the globally averaged radiative-equilibrium 
profile from Iro et al. (2005). Note the formation of a dayside inversion 
layer at pressures less than 0. 3 bars even though none exists in radiative 
equilibrium. 



ferences has increased from ~twofold (Fig. 4, left) to eightfold (Fig. 4, 
right). Emission from such high altitudes would largely mask the signa- 
ture of winds. It is possible that some planets have high- altitude clouds 
while others do not (HD209458b and TrES-1 may represent these cases; 
Fortney et al. 2005), so a range of infrared- lightcurve behaviors is to 
be expected among real planets. 

Dynamics can push the atmosphere far from radiative equilibrium, 
and this may have implications in explaining the existing Spitzer IRAC 
data at 4.5 and 8 /um for TrES-1 (Charbonncau et al. 2005). Current 
radiative-equilibrium models cannot easily explain the data: if they ex- 
plain the 4.5 /xm flux, then they do not predict enough flux at 8 /xm 
(Fortney et al. 2005, Burrows et al. 2005, Seager et al. 2005, Barman 
et al. 2005). Part of the problem is that 8-//m photons are emitted from 
higher altitude, where radiative-equilibrium models predict colder tem- 
peratures. One solution is to invoke a temperature inversion so that 
the emission region for 8-;um photons is hotter than for 4.5-/Ltm pho- 
tons. Fortney et al. (2005) accomplished this by adding an ad hoc heat 
source, which lead to an improved fit to the Spitzer data. A key point is 
that dynamics can naturally produce such a dayside temperature inver- 
sion, even when no such inversion would exist in radiative equilibrium 
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Figure 6.: Time evolution of the longitudinally averaged equatorial east- 
-west winds from the simulation in Fig. 1. Thin dashed, dash-dot, dot- 
ted, thick dashed, and solid curves show profiles at 1000, 2000, 3000, 
4000, and 5000 Earth days of simulated, time. Note the gradual down- 
ward penetration of winds from 10-100 bars over the course of the sim- 
ulation. 



(Fig. 5). The inversion occurs because of the upward-decreasing radia- 
tive time constant: as air columns superrotatc from nightsidc to dayside, 
the air at the top warms much more rapidly than air at the bottom, 
producing an inversion. No ad hoc heat sources need be invoked. 

The large radius of HD209458b has remained a puzzle, particularly 
given that the other 7 known transiting planets have radii in agreement 
with evolution calculations (e.g., Guillot 2005). Showman and Guillot 
(2002) and Guillot and Showman (2002) suggested that mechanical en- 
ergy produced by the atmospheric heat engine could be advcctcd into 
the interior, where it could be dissipated and might provide a source of 
heat that would slow the contraction. In the current simulations, all of 
the heating/cooling (which is the sole energy source in the simulations) 
occurs at pressures < 10 bars. Nevertheless, the simulations gradually 
develop strong winds at pressures > 10 bars, which implies that kinetic 
energy is transported downward from the heated regions into the inte- 
rior. Figure 6 shows the evolution of longitudinally averaged winds at 
the equator over time. The winds at pressures < Ibar rapidly reach 
a quasi-steady equilibrium, but the winds from 10-100 bars increase 
throughout the simulation (this implies a large increase in total kinetic 
energy because that layer contains ten times more mass than the entire 
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overlying atmosphere). The build-up of winds in Fig. 5 corresponds 
to a downward kinetic energy flux of lO^Wm-^, which is 10-100 
times greater than the intrinsic flux predicted in evolution models (e.g., 
Guillot and Showman 2002, Burrows ct al. 2003, Chabrier et al. 2004). 
More work is needed to determine the fate of this energy, but it suggests 
that atmospheric circulation could affect the long-term evolution. 
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